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We show that in supersymmetric models with a gauged flavor symmetry ex- 
plaining the hierarchy of fermion masses and mixings, the electron, muon, neutron 
and the deuteron acquire sizable electric dipole moments (EDM) through loop dia- 
grams involving the flavor gaugino/gauge boson near the Planck scale. These EDMs 
are proportional to the phases in the fermion Yukawa couplings and are typically 
much larger than the neutrino seesaw induced EDM for the leptons. In a popular 
class of models based on anomalous U(l) flavor symmetry of string origin, we find 
d e ~ (1(T 26 - 10 -27 ) e cm, df, ~ (10~ 24 - lO" 26 ) e cm, d n ~ 1(T 27 e cm, and 
dn ~ (10~ 26 — 10~ 27 ) e cm, which are within reach of next generation experiments. 



1 E-mail address: babu@okstate.edu 
2 E-mail address: enkhbat@okstate.edu 



1 Introduction 

The observed hierarchy in the masses and mixings of quarks and leptons is one of the 
most puzzling features of Nature. A plausible explanation is provided by flavor-dependent 
gauge symmetries. In a popular class of such models one extends the Standard Model 
(SM) to include a family-dependent U(l) factor. Upon spontaneous breaking of the 
U(l) symmetry effective Yukawa couplings of the form yij€ ni3 fif?H are induced, where 
e ~ 0.2 is a small parameter, and n^ are positive integers related to the family dependent 
U(l) charges. Even when the fundamental Yukawa couplings y^ are all of order one, a 
hierarchical spectrum is realized due to the suppression in powers of e Ij. Such flavor 
£7(1) symmetries can be naturally identified with the anomalous U(1)a symmetry of 
string theory [2]. Models using anomalous U(1)a symmetry for fermion mass and mixing 
hierarchy abound in the literature [2] 0. Most models of this type also assume low energy 
supersymmetry (SUSY) to stabilize the Higgs boson mass. Novel phenomena which are 
amenable to experimental tests can arise in such contexts. The purpose of this paper is 
to analyze one such effect, viz., the electric dipole moments of elementary fermions [7j. 

Low energy supersymmetry can potentially induce excessive flavor violation in pro- 
cesses such as K° — K° mixing and fi — > cy decay if the soft supersymmetry breaking 
Lagrangian takes its most general form. This potential problem is usually avoided by 
assuming a universal form for the soft SUSY breaking terms. Even with universality, the 
CP-violating phases present in the soft SUSY breaking Lagrangian can induce electric 
dipole moments (EDM) for the neutron and the electron at a level exceeding the current 
experimental limits. These effects have been extensively studied in the literature [%]-|14j. 
We will assume in the present work a universal SUSY breaking spectrum that is also 
CP-invariant so that excessive EDMs are not induced from the fundamental soft SUSY 
breaking parameters. 

The EDMs that we find in the context of models of fermion mass hierarchy are induced 
purely by complex Yukawa couplings. The phases in the Yukawa couplings are believed to 
be the source for the observed CP- violation in the K and B meson systems (CKM CP 
violation). It is thus reasonable to assume all Yukawa couplings, including the leptonic 
Yukawa couplings, to be complex. As we will see, it is natural that the flavor U(l) gauge 
symmetry responsible for explaining the fermion mass hierarchy breaks spontaneously at 
a scale Mp slightly below the fundamental Plank (or string) scale, Mp ~ M st /50. In the 
momentum regime Mp < fj, < M st the flavor gauge sector will be active and will contribute 
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to flavor violation in the squark and slepton sectors. These effects would survive down to 
the SUSY breaking scale and can lead to observable phenomena. In a previous paper we 
have studied leptonic rare decays \x — > cj and r — > /ry induced by the flavor gauge sector 
[E] . With complex Yukawa couplings, this flavor violation will also lead to EDMs for the 
electron (d e ), muon (d^), the neutron (d n ), and the deuteron (do) even with universal 
and CP-conserving soft SUSY breaking terms at the string scale. In a popular class 
of anomalous U(l) models which explains the fermion mass hierarchy, including bi-large 
neutrino mixing, we find d e ~ (10 -26 — 10~ 27 ) e cm, d^ ~ (1CT 24 — 10~ 26 ) e cm, d n ~ 1CT 27 e 
cm, and dp ~ (10~ 26 — 10 -27 ) e cm, which are within reach of next generation experiments. 
There are proposals to improve the current limit on electron EDM, \d e \ < 1.6 x 10~ 27 e 
cm [IS], by about two to four orders of magnitude [IBIHI]. It is expected that the current 
limit on the muon EDM, |d M | < 1.9 x 10~ 18 e cm, will be improved by six orders of 
magnitude or even more in the not too distant future JH] • There are also proposals which 
would improve the current neutron EDM limit from \d n \ < 6.3 x 10~ 26 e cm [TD] by a 
factor of 5 [2D]- The deuteron EDM is expected to be probed to the level of 10 -27 e cm in 
the near future |21j . Supersymmetry may reveal itself in these experiments before direct 
discovery at LHC, if the current ideas of solving the fermion mass hierarchy problem are 
correct. 

Lepton EDMs may arise even without flavor gauge symmetry from complex neutrino 
Yukawa couplings responsible for the seesaw mechanism in the context of low energy 
SUSY. This effect has received much attention recently [UJ 1221 122] • We have computed 
such effects for d e and d^, but found them to be much less significant compared to the 
flavor U(l) induced effects. For example, we find d e ~ 10~ 29 e cm for large tan/3 from the 
neutrino Yukawa coupling effects, to be compared with d e ~ 10 -26 e cm from the flavor 
£7(1) sector. Similar effects from GUT threshold has been studied in Ref. |24j . 

The paper is organized as follows. In Section 2 we review the class of models based 
on anomalous U(l) a symmetry for fermion mass hierarchy [H]. In Section 3 the EDMs 
induced by the flavor U(1)a gauge sector is analyzed. In 3.1 a qualitative discussion of 
the radiative corrections to the soft masses and A-terms is given. In 3.2 we present our 
full numerical results for the EDMs. Section 4 has our conclusions. In Appendix IA.1I we 
give the relevant expressions for the /3-functions for the soft SUSY breaking parameters 
including corrections from the U(1)a gauge sector. In Appendix IA.2I the fermion mass fit 
for the model used in the numerical analysis is presented. Appendix IA.3l lists the formulas 
needed for the calculation of EDMs. 



2 Fermion Masses and Anomalous U(l) Symmetry 



In this section we review briefly the idea of explaining fermion mass hierarchy with a 
flavor dependent U(l) symmetry. We focus on a specific class of anomalous U(1)a mod- 
els discussed in Ref. |H] to address the fermion EDM. Most models of Ref. [H] [2] will 
also fall into this category and will lead to similar results. In these models families are 
distinguished by their anomalous U(l) charges. The U(1)a symmetry is broken sponta- 
neously by an MSSM singlet flavon field S which acquires a vacuum expectation value 
(VEV) slightly below the string scale M st . This provides a small expansion parameter 
e = (S)/M st needed for explaining the fermion mass hierarchy. U(l) invariance forbids 
renormalizable Yukawa couplings for the light families, but would allow them through 
effective nonrenormalizable couplings suppressed by a factor (S/M st ) ni3 (for the fermion 
mass operator connecting flavors i and j) with n^ being positive integers. Even with all 
couplings being of order one, hierarchical masses for different flavors are naturally realized 
PP. Although this mechanism will work with any flavor U(l), anomalous U(l) models 
are attractive since they would also provide a natural understanding for the smallness of 
e ~ 0.2 [3], which arises from the one-loop induced Fayet-Illiopoulos D-term 
Consider the following fermion mass matrices studied in Ref jH]: 



M„ 



/«' 



M u ~ 



(H u ) 



e4 \ 



e 4 e 2 



e 4 e 2 



1/ 



M e ~ (H d )et 



fe 5 e 3 e\ 



1 
- 1/ 



e 4 e 2 



e 4 e 



(^ 



M 



R 



e 1 1 
e 1 1/ 



/e 5 e 4 



M d ~ (H d )e* 



I* 



M VD ~ 



(H u y 



M lig 



lit 



{ H u) _ e 2s 



e 1 
e 1 

2 



M, 



R 



^\ 



1/ 

£ \ 

1 , 
1/ 

1 1 
1 1/ 



(1) 



Here M u , M rf and M e are the up-quark, the down-quark and the charged lepton mass 
matrices in the basis fMff c (f = u, d, e, v). Complex order one coefficients multiplying 
each entry of the matrices are not shown. M VD and M u c are the neutrino Dirac and 
Majorana mass matrices. The light neutrino mass matrix M l J 9 is derived from the 
seesaw mechanism, p and s are integers and are chosen differently for different values 
of tan/? = (H u ) j (H d ). The choice of p = 0, 1, 2 corresponds to large (~ 50), medium 
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Table 1: The flavor U(1)a charge assignment for the MSSM fields and the flavon field S in the normal- 
ization where q s = — 1. In the third column we list the notation for the charges used in the paper. 



(~ 20) and small (~ 5) values of tan/3 respectively. The quark and lepton masses and 
mixings arising from Eq. (£Q) are fully consistent with experimental observations if e ~ 0.2. 
Note that the CKM mixing angles are small, while the leptonic mixing angles relevant for 
solar and atmospheric neutrino oscillations are of order unity. The lopsided nature of the 
matrices Md and M e of Eq. (0) enables such a disparity to be realized [2U] . 

The superpotential which would lead to the Yukawa coupling structure in Eq. ((TJ) has 
the following general form near the fundamental scale M st : 

+ fiH u H d + W A (S,X k ) . (2) 
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Here i, j = 1,2,3 are the generation indices and n{j (/ = u, d, e, v) are positive integers 
fixed by the choice of the U(1)a charge assignment. We choose all charges to be integers 
with the charge of S being negative. ?//• are the Yukawa couplings which we take to be 
complex and of order one for all i, j. H stands for the MSSM Higgs doublets H u and 
Hd- S is the MSSM singlet flavon field whose VEV determines the expansion parameter 
e. Wa in Eq. (J2J) contains MSSM singlet fields X k which would be needed for anomaly 
cancelation. The U(l) charge assignment shown in Tabled will lead to the texture of Eq. 
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The U(1)a anomalies are cancelled by the Green-Schwarz mechanism [2] which re- 
quires 
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Here A u A 2 , A 3 , A F and A gravity are U{lf Y x U{1) A , SU{2)\ x U{1) A , SU{3) 2 C x U{1) A , 



U(1) A and {Gravity) 2 x U(1)a anomaly coefficients. All other anomalies (such as U(1) A x 
U(1)y) must vanish, ki (i = 1, 2, 3), kp are the Kac- Moody levels, with the Non-Abelian 
levels k 2 and fc3 being integers. The factor 1/3 in front of the cubic anomaly Ap has a 
combinatorial origin owing to the three identical U(1)a gauge boson legs. 

We require string unification of all the gauge couplings including that of the U(1)a, 
gp, at the fundamental scale M st 



hgl = k F gp = 2g 2 st . (4) 

For a clear discussion of the coefficients in Eqs. ©-(0]) see Ref. J2H1- Here g^ are the 
U(l)y, SU(2)l and SU(3)c gauge couplings for i — 1, 2, 3. With the choice k 2 = k 3 = 1, 
consistency with the observed unification of gauge couplings within MSSM would require 
k\ = 5/3, corresponding to the 577(5) normalization of hypercharge. This will be an 
automatic consequence of the Green-Schwarz anomaly cancelation conditions with our 
choice of texture. The small discrepancy (in the absence of a covering GUT group) 
between the unification scale derived from low energy data and the string scale may 
be understood in the context of M-theory by making use of the radius of the eleventh 
dimension [2*9"] . 

With k 2 = h = 1 we find from Tabled A 2 = (19 + 3s)/2 and A 3 = (19 + 3p)/2. Eq. 
(JHJ) then requires p = s, i.e., a common exponent for the charged lepton and the neutrino 
Dirac Yukawa coupling matrices. With p = s, the condition A\/k\ = A 2 /k 2 fixes k\ to 
be 5/3, consistent with SU(5) unification. Note also that the charges given in Table d 
become compatible with S77(5) unification. Since Tr(Y) = for the fermion multiplets 
of S77(5), and since the Higgs doublets carry zero U(1)a charge, the anomaly coefficient 
[[/(l)^] 2 x U(1)y vanishes, as required. The last equality in Eq. (J3J) requires 

A gravity = Tr {q) = 12(19 + 3p). (5) 

This cannot be satisfied with the MSSM fields alone, since Tr(q)MSSM — 5(13 + 3p), 
which does not match Eq. (jSJ). We cancel this anomaly by introducing MSSM singlet 
fields X k obeying Tr (q) x = A gravity — Tr (q) MSSM = 163 + 21p. If all the X k fields have 
the same charge equal to +1, they will acquire masses of order M st e 2 through the coupling 
X k X k S 2 /M st and will decouple from low energy theory. For other choices of the charge 
of Xk these masses can be different. For example, if the charge is equal to +1/2, their 
masses will be of order M st e; if the charge is +2 the masses will be of order M st e 4 . We 
will consider only the case where the Xk fields have charge +1. 



With the charges of all fields fixed, we are now in a position to determine the U(1)a 
charge normalization so that g 2 F = g\ = g 2 at the string scale, (We take k 2 = k^ = 1.) 
This normalization factor, which we denote as |g s |, is given by \q s \ = 1/y/kp. All the 
charges given in Tableware to be multiplied by \q s \. From the Green-Schwarz anomaly 
cancelation condition Ap/^kp) = A2/&2, we have 

Tr (q 3 ) 19 + 3p 



2>kp 2k 2 

from which we find the normalization of the U(1)a charge |g s | = 1/y/kp to be 



(6) 



\q s \ = (0.179, 0.186, 0.181) for p = (0, 1, 2) . (7) 

The Fayet-Iliopoulos term for the anomalous U(1)a, generated through the gravita- 
tional anomaly, is given by [2H| 

t _ 9 s t M st \ i 4 / 8 \ 

^• — I no —2 I ys I ^gravity j \°J 

where g s t is the unified gauge coupling at the string scale (see Eq. (jlj)). By minimizing 
the potential from the U(1)a D— term 

v=^(^ r \sn^i\f a \ 2 +T,^\x k \A\ o) 

8 \\ls\ a k J 

in such a way that supersymmetry remains unbroken, one finds for the VEV of S 



e = (S)/M st = ^/g 2 st A gramty /192^. (10) 

For the fermion mass texture in Eq. (JT|). corresponding to the U(1)a charges given in 
Table HI we find 

e = (0.177, 0.191, 0.204) for p = (0, 1, 2) . (11) 

The masses of the U(1)a gauge boson and the corresponding gaugino are obtained from 
Mf = \qs\gF{S)/V2 and found to be 

M F =f^i,^,^forp=(0,l,2). (12) 

V54.5' 52.5' 53.9/ V ; V ; 

In the momentum range below M st and above Mp, these gauge particles will be active and 
will induce flavor dependent corrections to the sfermion soft masses and the A-terms. It is 
these effects which induce EDMs for the electron, muon and the neutron at low energies. 
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3 Electric Dipole Moments from Anomalous U(l) 

In the Standard Model the electric dipole moments of the electron, muon and the neutron 
are predicted to be extremely small and beyond reach of planned experiments. In the 
presence of low energy supersymmetry these EDMs can exceed the current experimental 
limits if soft SUSY breaking parameters are complex [H] [T3|. To focus on the anomalous 
U(l) induced effects we shall adopt the minimal supergravity scenario with universal 
and CP-conserving soft SUSY breaking parameters. Specifically, at the string scale we 
assume a universal scalar mass m , a common gaugino mass M\/2, and trilinear A-terms 
proportional to their respective Yukawa couplings. We assume m , M\ii and A , and 
the Higgs mass parameters ft and Bfi to be real. Thus the only source of CP-violation 
is in the complex Yukawa couplings. This is needed for the CKM CP-violation in the 
quark sector and it is natural to assume that the leptonic Yukawa couplings are complex 
as well. In Sec. 3.1 we give a qualitative estimate of the EDMs induced by the radiative 
corrections involving the U(1)a gauge sector in such a SUSY context. Our numerical 
results are presented in Sec. 3.2. 

3.1 U(1)a Correction to the Soft Parameters and EDM 

We now give approximate expressions for the U{1)a gauge sector RGE corrections to the 
soft parameters between the string scale and the U(1)a breaking scale Mp. The full RGE 
expressions for the soft parameters in the presence of higher dimensional operators as 
in Eq. (J2J) have been derived in Ref . JO] . In Appendix IA.1I we summarize the relevant 
expressions. The U(1)a corrections to the soft masses for the left-handed slepton are 
obtained from Eq. (|2*T|) to be 

5 (m?);* * (4(gf M Xp f - «Tr (<*)) (ksM% l ° g (M J Mp) , (13) 

and a similar expression for the right-handed slepton masses with the interchange (L, q L ) 
— > (e, q e ). There are analogous corrections in the squark sector. The corrections to the 
A-terms are obtained from Eq. (J2"^j) as 



oAij ~ -Mx F g F Yij^ij 4-5 , ( 14 ) 

where Z^ are biliear combinations of the flavor charges given by [0] 

Ztj = qN] + qrh + qfh + ny s (q^ + q e j +h) + -n'i j (n^-l)q 2 s . (15) 



Numerical values of Zf- for our model are given in Eq. (J26J1 of Appendix IA.lt Note that 
these corrections, Eqs. (fT3"j) and (JT3J), are flavor dependent. Due to the flavor dependent 
nature of these corrections, the fermion and the corresponding sfermion mass matrices 
cannot be diagonalized simultaneously. This was the source of the flavor violation studied 
in Ref. j^j. For the same reason, with complex Yukawa couplings Y/j, nonzero EDMs for 
the fermions will be induced. 

Let us now estimate the EDM of the electron arising from the corrections in Eqs. 
()13j) and p4jl . There are three flavor dependent matrices in the leptonic sector, not 
including the neutrino Yukawa matrix Y u . They are the leptonic Yukawa matrix Y e and 
the matrices of U(1)a charges Q L = diag (1 + p, p, p) and Q e = diag (4, 2, 0) for the 
lepton doublets and singlets (see Table HJ). In the mass eigenbasis for the charged leptons 
Q L and Q e will develop complex off diagonal entries, with the phases arising from Y e 
through the unitary matrices that diagonalize Y e . This is the basic source for the EDM. 
The corrections given in Eq. (|13|) will generate EDM of the electron through the product 
of slepton mixings in (H)ll, (H)lr and (H)rr (for % = 2,3). The induced EDM will 
be d e ex Im (U^Q L Y e Q e V H , where U and V are unitary matrices which diagonalize 
Y e , Y e = UYj ai V'. There are additional corrections which are quadratic in Q L and Q e . 
The corrections to the A-terms in Eq. (|T4"|) will also induce EDM directly through (LR) 
mixings. Combining these effects with the formula for the EDM given in Eq. (|42j) of 
Appendix IA.3I we arrive at the following approximate expression for d e : 

~V e V e 



' 8vr cos 2 9 W m\ \m\) 8vr 2 ^ L l 



Im 
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it 



(16) 



where C™ and Cf denote the contributions from the soft masses and the A-terms respec- 
tively. They are given by 

_ (\q s \g F ) 2 \og(M st /M F )m 4 (A - Htan/3) ^ Lrji? 

iff = 4 (M 1/2 /m ) 2 ((gf f - (gf) 2 ) - (gf - gf)Trg, 

Cf = 2^ {Zt, - Z{ x ) . (17) 

m i 

Here iff is obtained from iff by the replacement gf — > g| . mj is the average slepton 
mass and M^ is the Bino mass. The function A (X) is given in Eq. (|43|) in Appendix 
IA.31 We see explicitly that the complex Yukawa couplings along with nonuniversal U(1)a 
charges lead to nonzero EDM. 



To estimate the size of this effect we choose the approximations m = M 1 / 2 — Msusy- 
Following the mass matrices given in Eq. (£Q) we take \Y2\ ~ e n ^ +p . We consider here 
only the contribution from the (13) mixing, since the U(1)a charge difference is the largest 
between the first and the third generations in Q e . Then we find 

a / An-27 \ /500GeV\ 2 

d e /e ~ (10 ^'cm x 



x mJo(io) ^^P^ + o { i)M^l) Alg[yiWll . (18) 

\ m l m l I 

From this estimate we see that the electron EDM induced by the U(1)a gauge corrections 
is in the experimentally interesting range and already puts constraint on the soft SUSY 
breaking parameters. The actual numerical result is quite sensitive to the choice of ttiq 
and Mi/2- In our numerical calculations we have chosen m = Mi/2/4.4 for low tan/3 for 
cosmological reason. In this case the O(10) coefficient in Eq. (fTK|) will be reduced to an 
0(1) number. For large tan/5 this coefficient will remain as 0(10). 

Let us now compare the anomalous U(l) induced EDM with the right-handed neutrino 
induced effects pointed out in Ref. J3] and studied further in Ref. J22I23- The latter 
effects induce EDM which are given by 

dt ocKrf Ay, (rfY"], (19) 

where Ay = \og{McuT I '{M v c)i)5ij . Here (M„ c )j is the mass of the right-handed neutrino 
of flavor i. With our texture for the neutrino mass matrices dictated by U(1)a symmetry 
we find the right-handed neutrino induced EDM to be d e ~ 10~ 29 e cm, which is two 
to three orders magnitude smaller than the anomalous U(1)a induced effects. In our 
numerical analysis we present separately our results for the electron EDM arising from 
the right-handed neutrino effects. 

3.2 Numerical Results 

In this sub-section we present our numerical results for the electron, muon, neutron and 
the deuteron electric dipole moments. We adopt the minimal supergravity scenario for 
supersymmetry breaking. At the string scale, taken to be M st = 10 17 GeV, we assume 
a universal scalar mass mo and a common gaugino soft mass M1/2. All SUSY breaking 
parameters (mo, Min, A , B ) and the \x term are taken to be real at the string scale. 
We choose /x > for all cases except in Fig. ^ where we also show results for /i < 0. 
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The anomalous U(l) gauge coupling g F is chosen to be gp/4:7i = 1/24, consistent with 
string unification. The soft SUSY breaking parameters are evolved from M st to the U(1)a 
breaking scale Mp ~ M st /50 (see Eq. (fl2"jl ) including the U(1)a gaugino/gauge boson 
corrections. 

We present our results for the EDM for three values of the parameter tan/3, small (5), 
medium (20) and large (50). We take m = Mj/2/4.4 for low and medium values of tan/3. 
This is motivated by the requirement that the right abundance of cosmological cold dark 
matter be generated. With mo = Mi/2/4.4, the right-handed charged sleptons will have 
masses slightly above that of the neutralino LSP. The relic abundance of neutralinos is in 
the right range with such a spectrum, as a result of coannihilation [3U|. For large tan/3 we 
also allow the choice mo = M1/2, since alternative mechanisms for reproducing the right 
relic abundance of LSP become available in this case [TTTj . 

We vary Mi/ 2 in the range 250 GeV to 1 TeV. The results are presented for two 
different values of A , and 300 GeV. The lepton EDMs induced by the flavor U(l) 
gaugino/gauge boson contribution are plotted against the universal gaugino mass M1/2 in 
Figures QE1 In Figure El the electron EDM induced purely by the right-handed neutrino 
effects is plotted. In Figure EHH1 we plot the EDM of the neutron and the deuteron arising 
from the flavor U(l) gauge boson/gaugino effects. 

As input at M st we choose the Yukawa coupling matrices given in Eqs. (|37|) - (|4lH) 
of Appendix IA.2I (for tan/3 = 5). These are obtained by extrapolating the low energy 
Yukawa couplings to M st and applying bi-unitary transformations at M st to generate the 
texture given in Eq. ( jlj) . The low energy Yukawa couplings and their extrapolation are 
discussed in Appendix IA. 21 As for the neutrino Dirac Yukawa couplings, we choose Y v to 
be such that in the flavor basis (after the bi-unitary rotations) it exhibits approximately 
the structure given in Eq. ([!} with (Y u ) 33 ~ e p . For a given choice of hierarchical 
light-neutrino spectrum this would uniquely fix the right-handed neutrino mass matrix 
through the seesaw mechanism. M v c will then have the form given in Eq. (|I|) . We set 
(M^c) 33 = Mft6 2p with M^. ~ 4 x 10 14 GeV. The eigenvalues of the right-handed neutrino 
mass matrix are important for the lepton EDMs induced by the right-handed neutrino 
threshold effects. It should be noted that the unitary rotations applied on the diagonal 
Yukawa matrices at M st are not unique, except that they should conform to the fermion 
mass matrix structure shown in Eq. ([TJ. So our fits should be taken only as indicative, 
and not definitive. We expect differences of order one in our numerical results on EDM 
arising from the arbitrariness in these unitary matrices. 
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In Figured] the electron EDM induced by the U(1)a gaugino/gauge boson contribu- 
tions to the soft masses and A-terms are plotted as a function of M1/2 for three values of 
tan /3. We see that some parts of the parameter space are already excluded by the current 
experimental upper bound d e < 1.6 x 10~ 27 e cm and that the other parts are in the range 
which will be tested by next generation electron EDM experiments |16j . 

In Figure 121 we plot the muon EDM as a function of SUSY breaking parameters. We 
find d^ to be in the range (10~ 25 — 10~ 28 ) e cm for most of the parameter space. This 
value is somewhat smaller than than d e (m M /m e ), which would be the naive expectation 
based on the scaling of lepton masses. This happens for the following reason. The second 
and third family left-handed charged sleptons have the same U(1)a charge, so the flavor 
gauge bosons/gaugino will not generate any mass splitting between these sleptons. The 
mixing in the right-handed charged slepton sector is suppressed by a factor e 2 for all tan f3, 
compared to the suppression factor e between the first and the second generations. On the 
other hand, we find quite an enhancement of the muon EDM for the choice m = M\i<i and 
tan/? = 50. For this choice, the electron EDM is well above the experimental bound. Since 
the two EDMs are induced by independent phases, it is possible to choose the parameters 
such that the electron EDM is below the experimental limit and at the same time the 
muon EDM is at the level of ~ (1CT 25 — 10 -24 ) e cm, although we do not attempt such 
an explicit solution here. It should also be pointed out that parts of the parameter space 
where d^ is large is already ruled out by the experimental upper limit for the radiative 
decay fi — > c-f for the numerical fits shown jH]. The remaining regions will be put to 
experimental scrutiny by future experiments [T8] . 

In Figure El we present for comparison, the electron EDM arising solely from the right- 
handed neutrino threshold effects [14 . With the proper decoupling of the right-handed 
neutrinos [22] we find our results to be in rough agreement with those in Ref. [T3| 1221 123! - 
Nevertheless these effects, which yield at most d e ~ 10~ 29 e cm, are much smaller compared 
to the U{1)a effects. 

In Figure El we plot the neutron EDM versus M1/2. In Figure El we plot the deuteron 
EDM. Details of the calculations are given in Appendix IA. 31 In both cases our numerical 
results are in the interesting range which should be accessible to proposed experiments in 
the near future. We find the contributions from the CKM phase to be of the same order as 
the contributions from the U(1)a gaugino/gauge boson sector. Figures EHI3 include both 
these effects. The flavor sector contribution to the neutron EDM is somewhat smaller 
compared to the leptonic EDM due to the gluino focusing effect. (The squarks receive 
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flavor universal contributions for their masses below Mp from the gluino, which tends to 
suppress flavor violation and thus d n .) 

We have also studied the constraint on the chromoelectric dipole moment for the 
strange quark dg arising from 199 Hg EDM [321 ES]. This bound reads as \d^\ < 5.8 x 
10 -25 e cm. This constraint is easily satisfied in our model. The down-type squark mixing 
in the (23) sector is suppressed by a factor e 2 for the right-handed squarks, and is van- 
ishing to leading order for the left-handed squarks, similar to the case of \x — r mixing. 
Consequently, we find the chromoelectric EDM of the strange quark to be about two to 
three orders of magnitude below the experimental limit. 

The soft SUSY breaking bilinear 5-term and the gaugino masses will develop complex 
phases via the one-loop and two-loop RGE corrections respectively arising from the A— 
term contributions. In our model we find these corrections to be negligible compared to 
the U(1)a flavor gaugino/gauge boson effects. 
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Figure 1: Electric Dipole Moment of the electron induced by the flavor gaugino/gauge boson. The 
(red) horizontal line shows the current experimental limit on d e . We have chosen here too = M 1 / 2 /4:A. 
For tan/3 = 50 we show an additional case with Too = M 1 / 2 (the uppermost curve). For tan/3 = 20 and 
Ao = 300 GeV we find a cancelation between the A-term contributions given in Eq. (|24|l and the soft 
left/right mass contributions in Eq. (|21|) for our particular fit of the Yukawa couplings. This cancelation 
disappears for the choice of negative /z-term (the curve labeled by /x < 0) . 
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4 Conclusions 

In this paper we have studied the electric dipole moments of the electron, muon and the 
neutron induced by a flavor dependent U(l) symmetry which explains the hierarchy of 
fermion masses and mixings in a natural way via the Froggatt-Nielsen mechanism. This 
£7(1) symmetry may be identified as the anomalous U(l) of string theory. This symmetry 
is broken spontaneously at a scale Mp slightly below the string scale, Mp ~ M st /50. In 
the momentum regime Mp < \x < M st , the flavor U(1)a gauge boson sector will be active 
and will contribute to the soft SUSY breaking parameters in a flavor dependent fashion. 
We adopt the minimal supergravity scenario for SUSY breaking, and assume that the 
soft SUSY breaking parameters are universal and real. The complex Yukawa couplings 
will still induce phases in the soft SUSY masses and the A-parameters, leading to the 
generation of EDM. This is the main source of the EDM that we have studied here. 

We have presented our numerical results for the electron, muon, neutron and the 
deuteron EDMs in Figures EJH1 as functions of supersymmetry breaking parameters. d e 
and d n are very close to the current experimental limits, d e ~ (10 -26 — 10~ 27 ) e cm 
and d n ~ 10~ 27 e cm. For the deuteron, our prediction, dp, ~ (1CT 26 — 10~ 27 ) e cm, 
would make it within reach of proposed experiments. For the case of the muon, although 
d^ is rather small for low tan/3, in the case of large tan/3 ~ 50, for certain choices of 
phases in the Yuakawa couplings, we have found the induced EDM to be as large as 
d^ ~ (10~ 25 — 10 -24 ) e cm, which might be accessible to future experiments [IS] . In 
the leptonic sector, these EDMs are much larger than the ones induced by the neutrino 
seesaw sector, which yields, for example, d e ~ 3 x 1CU 29 e cm with our texture of fermion 
mass matrices dictated by flavor symmetries. In Figure H] we present our results for the 
induced d e arising from the neutrino seesaw sector. Discovery of electric dipole moments 
for the electron, muon and the neutron can shed light on one of the fundamental puzzles 
of Nature, viz., the origin of mass for elementary particles. 
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A Appendix 

A.l RGE for Soft Parameters including U(1)a Corrections 

Here we list all the /3-functions for the soft masses and A-terms including contributions 
from the anomalous U(l) gaugino/gauge boson sector (see Ref. for details). 
The one-loop /3-functions for the soft sfermion masses are given by 

MSSM / \A / \N 

) +P(m 2 f ) +P[m 2 f) 

' ij \ ■) / ij \ J / ij \ J / ij 

where the superscripts MSSM, A and N denote the contributions from the MSSM sec- 
tor, the anomalous U(l) gauge sector, and the right-handed neutrino sector respectively. 
/3 I m f are present only for the left-handed slepton and the sneutrino soft masses and 

V J/ij 

their explicit form is given elsewhere. The anomalous U(l) part (5 (m|) is [35. 

V J/ij 

/3 (rrff = ^-^iglh (cr - 4g/(M AF ) 2 ) . (21) 



/ o\ / n \MSSM / „\A / „\JV , 

0(m%=0 H)„ + I3{mj) + fj(mj), (20) 



ij 167T 

Here a is defined as 



a = 3 Tr \2q^m z Q + q u m z u + q a rhj) + Tr \2q L m l L + q e m z e + q v m{ 

+q s m 2 + J2Q?rh 2 x h > ( 22 ) 

k 

where fhx k are the soft masses of the extra particles X^ introduced for anomaly cance- 
lation via the Green-Schwarz mechanism. The trace is taken over family space. Here 
P(m 2 L ) i : ' etc., stand for the MSSM /5-functions without the neutrino and the flavor 
U(1)a contributions. 

Introducing a notation 

4 = 4^. ( 23 ) 

the U(l) gaugino part of A-term /3-functions is given by 

Pi£*)* = ~^9 2 f4 {(q!) 2 + (if) 2 + h 2 ) + ^g 2 F Z{ 3 Y^M Xf , (24) 

where g F and M\ F are the U(1)a gauge coupling and the gaugino mass respectively, h is 
the U(1)a charge of the up-type (down-type) higgs doublet if f c is up-type (down-type). 
Here we defined the combination of the U(1)a charges Z(, as 

Z lj = q{qf + qlh + qfh + n{ j q s (q{ + qf+h) + -n{ j (n{ j -l)q 2 . (25) 
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From the U(1)a charge assignments for the MSSM fields as given in Table [l"]one has 

/(ll, 13, 16) (4,6,9) (1, 3, 6)\ 



Z e = - 



(10,11,13) (3,4,6) (0,1,3) 
V(10, 11, 13) (3,4,6) (0,1,3)/ 



/(ll, 13, 16) (9,10,12) (9, 10, 12) \ 

(4,6,9) (3,4,6) (3,4,6) 

V (1,3,6) (0,1,3) (0,1,3) / 



(26) 



(27) 



/(2, 4, 7) (1,3,6) (1, 3, 6)\ 
Z v = - (1,2,4) (0,1,3) (0,1,3) , (28) 

V(l, 3, 4) (0,1,3) (0,1,3)/ 

for the three different values of p= (0, 1, 2) and 

/20 13 10\ 

Z u = - 13 6 3 , (29) 

\10 3 / 

independent of p. 

The terms we are interested in are the ones proportional to the U(1)a gaugino mass 
M\ F and the term proportional to a in the /3-functions. Beside these, the A-term j3- 
functions contain a flavor dependent piece which arises from the wave-function renormal- 
ization. Since the corresponding Yukawa /3-functions contain the same terms, these are 
simultaneously diagonalized, and do not lead to flavor violation. 



A. 2 Fermion Mass Fit 

Here we present the numerical fits to the fermion masses and mixings adopted for the 
calculation of the EDMs. As input at low energy we choose the following values for the 
running quark masses [3*4*] : 

m„(lGeV) = 5.1 MeV, m c (m c ) = 1.27 GeV, m t (m t ) = 167 GeV, 

m d {l GeV) = 8.9 MeV, m s (l GeV) = 130 MeV, m b (m b ) = 4.25 GeV. (30) 

The CKM mixing matrix is chosen in the standard parametrization with # 12 = 0.221, # 13 = 
0.005, #23 = 0.043 and the complex phase S = 0.86. We Use two-loop QED and QCD 
renormalization group equations to evolve these masses from the low energy scale to the 
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scale of SUSY breaking. We obtain the following running factor rf = m f (Msusy) /'n^fi'm'f) 
for the fermion masses at the SUSY breaking scale MguSY, initially chosen to be 500 GeV 
with a s {M z ) = 0.118: 



{r t , r b , r T , r u , r c , r d , s , r e , M ) = (0.943,0.605,0.991,0.395,0.442,0.395,0.398) 



(31) 



Then these masses at Msusy are used to calculate the Yukawa couplings in DR scheme. 
Using one-loop SUSY RGE evolution above Msusy we obtain the Yukawa couplings at 
the U(l) A breaking scale (M F ~ 10 15 GeV) to be 

(Y u , Y c , Y t ) = (5.2699 x 10~ 6 , 1.4634 x 10" 3 , 0.55498), 

(Y d , Y s , Y b ) = (3.4415 x 10" 5 , 5.0222 x 10~ 4 , 2.8247 x 10~ 2 ), 

{Y e , Y„, Y T ) = (1.0388 x 10~ 5 , 2.1481 x 10~ 3 , 3.6239 x 10" 2 ), 

(Y Ul , Y V2 , Y V3 ) = (7.4112 x 10~ 4 , 3.5135 x 10" 3 , 4.7212 x 10~ 2 ), (32) 



for tan/? = 5, 



for tan /? = 20 and 



(Y u , Y c , Y t ) = (5.0919 x 10~ 6 , 1.4140 x 10~ 3 , 0.53090), 
(Y d , Y s , Y b ) = (1.3834 x 10" 4 , 2.0189 x 10~ 3 , 0.11508), 
(F e , Y^ Y T ) = (4.1778 x 10~ 5 , 8.6456 x 10~ 3 , 0.14818), 
(Y Ul , Y U2 , Y U3 ) = (3.5720 x 10~ 3 , 1.6934 x 10~ 2 , 0.2275), 

(Y u , Y c , Y t ) = (5.3161 x 10~ 6 , 1.4764 x 10~ 3 , 0.59610), 
{Y d , Y s , Y b ) = (4.2226 x 10~ 4 , 6.1621 x 10" 3 , 0.41186), 
(Y e , Yp, Y T ) = (1.2720 x 10" 4 , 2.6645 x 10~ 2 , 0.51807), 
(Y Ul , Y U2 , Y U3 ) = (1.7822 x 10~ 2 , 8.4492 x 10~ 2 , 1.1354), 



(33) 



(34) 



for tan/3 = 50. 

We have chosen the Dirac neutrino Yukawa couplings such that in the flavor basis 
(after the bi-unitary rotations) it exhibits the assumed structure in Eq. (£Q). The light- 
neutrino masses and the lepton mixing matrix are chosen to be m vi = 2.7 x 10 -3 eV, 
m Vr , = 6.4 x 10" 3 eV and m„, = 8.6 x 10~ 2 eV and 

/ 0.8494 -0.5262 -0.04 \ 

0.3915 0.5775 0.7164 , (35) 

V -0.3539 -0.6242 0.6965/ 



'-(/_> 



''V-A 



V, 



MNS 



which are compatible with the observed light-neutrino oscillation parameters. From this 
the right-handed neutrino mass matrix is fixed uniquely with (M u c) 33 ~ M^e 2p where 
Mft = 4 x 10 14 GeV. The eigenvalues of M v c are used to calculate the lepton EDM 
induced by the right-handed neutrino threshold effects. 

In the following, we present our fits to the texture of Eq. ((TJ which have been used 
in our numerical calculations for tan/5 = 5 (We have similar fits upto an overall factor 
for tan/3 = 20, 50). This fit is not unique, and so can lead to order one uncertainty in 
our EDM results. The following fit is found by applying bi-unitary transformations with 
complex phases on the diagonal Yukawa coupling matrices at the U(1)a breaking scale 
Mp. We introduce a notation for the Yukawa couplings: 



-/ 



YA = y^e-v 



it A 



(36) 



At Mp we find the following fit (for tan/? = 5): 

( (1.45 + 1.6(H) e 8 (-0.563 -1.24i)e 6 

;-0.769 - 0.584i) e 6 (0.765 - 0.109i) e 4 



Y 



V (-0.282 - 0.204i) e 4 (0.274 - 4.40 x 10~ 2 i) e 2 



(1.50-0.397i)e 4 

;-0.255- 0.261 x 10~ 2 i)e 2 
0.554-2.80 x 10" 5 i 



(37) 



Y" 



( 



1.87- 1.69 i)e 5 



1.93 + 0.849 i)e 4 



1.29 + 0.957i)e 4 \ 



-0.404- 0.248 i)e 3 (0.552 + 1.54 x 10~ 2 z) e 2 (0.702 - 0.546 i) e 2 
V (-1.52 -0.435 i)e 0.312 + 0.314i 0.543 - 4.74 x 10~ 4 i/ 



(38) 



Y' 



/(3.52x 10- 2 + 0.480 i)e 5 

(-0.170 -0.612i)e 4 
V (0.538- 0.421 i)e 4 



(-1.85- 1.74 i)e 3 

1.15- 4.65 x 10- 2 i)e 



-0.539- 0.579 i) e 
0.319- 0.321 i 



-0.419 - 0.536 i) e 2 0.784 - 9.73 x 10~ 4 i J 



(39) 



/ (0.232- 0.190 i)t 2 (0.217-6.09x10 



-2 



i e 



Y' 



-0.206 -0.637i)e \ 
(0.638- 0.652 i) e -7.82 x 10" 2 + 0.537i 0.804 + 0.296 i 

\ (0.305 -0.392 i)e -4.41 x 10" 3 + 0.277 i 0.404 - 3.89 x 10 2 i / 

Note that all coefficients multiplying e niJ in Eqs. fj3T|) - (|4T?|) are of order unity. 



(40) 



A. 3 Formulas for Electric Dipole Moments 

We list here the formulas for the electric dipole moments of leptons and quarks in the 
MSSM from Ref. JUJj which we have used in our numerical analysis. 
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The EDMs of elementary fermions are sum of neutralino, chargino and for quarks 
gluino contributions which we denote as d^, dfj and d^. In addition to these, the quarks 
receive contributions from chromoelectric and purely gluonic dimension-six operators [36] . 
We have not considered here the latter one, since these effects turn out to be small. The 
effective EDM operator df for a spin-| particle is given by 



C = --dfiPa^F^ 



(41) 



The EDM df in general has the following components in a supersymmetric theory: 

m , 



d%/e 



a 



8n sin 2 9 



w x=1 a=1 m f ~ x 



d C Je 
d c Je 



-a 



6 



87r sin 9 W 
—a 



EEMQ, 



M, 



*h 



x=l b=\ 
6 2 



b 2 

d x 



8vr sin 2 9 



/ , \ M v+ / 2 



W x =l b=l 



ml_ 




«.s 



<■■> 



M~ l M -+' 

-^EM^r '"-"' Xb 



3n 



x=l 



m% 



-Q~ qi A 



2 



(42) 



where 



A(X) 
B(X) 



l-X 2 + 2XlogX 

(1 " ^) 3 ' 

3-4X + X 2 + 21ogX 



(43) 



Here M^o, M-+ and M^ are the neutralino, chargino and the gluino masses respectively. 
m 2 ; (x = 1, ..., 6) are the eigenvalues of the sfermion mass matrices. The coefficients iV"/ , 
C xh and G* are given by 

Nil = {V2tan9 w Q fi (0 N ) ia U! +3iX -K f (0 N ) a)a Ul x 

x -V2{t a n9 w (Q fi - T h ) (o*)^ + T 3/i (o") J flg 






tf«(og)>L 



62 



o 



F) 61 ^t-^(o£) 62 C/t 



^K) fc >-[K) W ^-^K), 2 ^ 



7-r? 779* 



(44) 



where -K" u = m u / (\^2M W sin (3) and Ki t d = mi >( i/(V2M w cos/5). O^ and 0£, O^ ma- 
trices diagonalize the neutralino and chargino mass matrices respectively. The index a' 
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of N in the neutralino contribution formula takes value of 3(4) for T 3 j = — |(|). The 
chromoelectric dipole moments d q for quarks are defined as 

Ccedm = - % -g s d q qT a a, vlb qG^ a . (45) 

The contributions to d q from neutralino, chargino and gluino are given by 



where 



C{X) 



Zl^tih ™l [ml 

_„2 6 2 M-+ (M 2 . + 

da = ^EE Im (C| 6 )HM 4 

a 6 M- ( m2 +\ 

€ = ^EMG^^MM' (46) 



f 10 X-26 + ^£_^l£V (47) 



6(1-X) 2 V^ ~" ' 1--^ !"^ 

We use the QCD sum rule based estimate of Ref. [37] to evaluate the neutron and the 
deuteron EDMs: 



d n = 0.7 (d d - 0.25d u ) + 0.55 (d d + 0.54) , 

d D = 0.5 (4 + 4) - 0.6 (4 - 4 + 0.3 (d d + 4) ) . (4* 

Here the running factors are d q (1 GeV) ~ 0.914 (M z ) and d q (l GeV) ~ l.2d q (M z ). 
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